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Abstract 

A Quark-Meson Coupling (QMC) model is extended to finite nuclei in the 

relativistic mean-field or Hartree approximation. The ultra-relativistic quarks 

are assumed to be bound in non-overlapping nucleon bags, and the interaction 

between nucleons arises from a coupling of vector and scalar meson fields 

to the quarks. We develop a perturbative scheme for treating the spatial 

nonuniformity of the meson fields over the volume of the nucleon as well as 

the nucleus. Results of calculations for spherical nuclei are given, based on a 

fit to the equilibrium properties of nuclear matter. Several possible extensions 

of the model are also considered. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD), a theory of quarks and gluons, is the current 
paradigm for the strong interaction. Nuclei are bound by the strong interaction, but hadronic 
degrees of freedom account for the most of the physical properties. This contrast suggests 
both questions and opportunities. One might ask how and why this standard picture of 
nuclear physics as a system of clustered color-singlet objects emerges from the fundamental 
theory. The opportunities arise from the challenge of discovering small (at normal nuclear 
densities) but interesting corrections to the standard picture which arise from the underlying 
constituents. One example was the discovery of the EMC effect [1] , in which the structure 
function of a nucleon was shown to be modified by the nuclear medium. 

Thus the need for a theory of nuclei that incorporates quark-gluon degrees of freedom, but 
also respects the vast body of information substantiating the standard picture, is apparent. 
There are many possibilities, and here we wish to focus on the Quark-Meson Coupling 
(QMC) model. 

In a recent series of papers Saito and Thomas [2] have considered a model for nuclear 
matter involving color-singlet clusters of quarks and also mesons. This model is a variation 
on one originally proposed by Guichon [3], with corrections by Fleck et al. [4]. The quarks 
are assumed to be bound in non-overlapping nucleon bags, and the interaction between 
nucleons arises from a coupling of meson fields to the quarks. The ultra-relativistic quarks 
are described by a mean-field Dirac equation together with the MIT bag model boundary 
conditions [5]. The nucleons are assumed to be also described by the Dirac equation in the 
effective mean fields arising from meson fields coupling to the quarks in the nucleon. The 
model is conceptually similar to the QHD model of Walecka and collaborators [6] , in which 
the meson fields couple to point-like nucleons. Indeed, in the limit that the quark mass 
becomes very large the present model reduces to the Walecka model. But the presence of 
the quark degrees of freedom provides important physical content absent from QHD. 

The purpose of the present paper is to outline the formalism for applying this model to 
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finite nuclei. Justifications for the adiabatic approximation we use, in which the quarks do 
not respond to the motion of the nucleon within the nucleus, have already been obtained 
by Guichon et al. [7]. These authors also reported preliminary results for ^^O within a 
local density approximation to the nuclear matter results. Here we provide a more detailed 
formalism which avoids the need for some of the approximations they use. We also discuss 
extensions of the model to account for other effects that may depend on the nuclear density 
[8] , and the use of a relativistic oscillator description of the nucleon [9] as an alternative to 
the bag model. 

An outline of the paper is as follows. The next section is concerned with the formalism 
for both finite nuclei and nuclear matter. Results are then given for spherical nuclei, based 
on a fit of the model parameters to the properties of nuclear matter. Extensions of the 
model are then considered, followed by a summary and outlook in the conclusions. 



II. THEORY OF FINITE NUCLEI 
A. Nuclear Energy and General Field Equations 

The nucleus is treated as a set of color singlet objects carrying nucleonic quantum num- 
bers. These objects (which we shall denote as bound nucleons) are not the same as free 
nucleons because the forces that bind the nucleus also infiuence the baryonic wave function. 

In the present model the forces between the bound nucleons are generated by the ex- 
change of vector (a;) and scalar (cr) mesons coupling to quarks. The static Hartree approx- 
imation is used to describe the meson mean-fields and the bound nucleon single particle 
states. 

As in the QHD model [6], we will assume the nucleons obey the Dirac equation 

" -z7-V + M*(X)+70K(X)]^,(X)=S,7ViW, (1) 

with nucleon states labeled by quantum numbers i. The quantity M*[X) is the mass of the 
bound nucleon, which differs from the mass of the free nucleon due to the infiuence of the 



scalar potential generated by the sigma field on the quarks inside the nucleon. As will be 
shown, M*{X) is a nonlinear function of the sigma field in the QMC models. Vy{X) is the 
time component of the vector potential, and depends linearly on the time component of the 
omega field. The vector component of the omega field averages to zero in the nuclear rest 
frame, so we drop it here. 

The total nuclear energy S is obtained as a sum of the nucleon and meson field energies: 

£^ Ei + ^v + Es. (2) 

i=occ 

The first term is a sum over the occupied positive energy nucleon states. The contribution 
of the vector mesons to the total energy is given by 

Ev = -\j d'r a;(r) {ml - V^) a;(r), (3) 

and the sigma meson contribution is 

Es^\ j d\ a{r){rnl-V'')a{r). (4) 

The mass of the vector meson is taken here to be the same as the omega (m^=783 MeV), 
whereas the mass of the scalar meson is a parameter of the model, typically in the range 
rus ~ 400 - 550 MeV. 

The fields ijj{r) and a{r) provide vector 

V,^{r)^glu^{r), (5) 

and scalar quark potentials 

VKr)=gl<j{rl (6) 

that act on confined quarks located at r. Here and gl are the relevant quark-meson 
coupling constants. The quark potentials determine the quark wave functions and densities, 
which in turn are needed to obtain the functions oj{r) and (T{r). Thus there is an additional 
self-consistency requirement in the present model — the quark and nuclear wave functions 
must be determined together. 
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The meson field equations are derived by minimizing E with respect to variations in the 
meson fields a;(r) and cr(r), yielding 

(-V^ + m^) u^ir) = / d^X P.{X)^-j^, (7) 
(-V^ + m^) a{r) - - J d'X Ps{X)^-^^. (8) 

Here p^^siX) represents the vector (v) and scalar (s) nucleon densities at a position X, with 



p,{X) = J2 ^^{X)^P^{X), (9) 

i=occ 

Ps{x) = y: 4'^{xmx). (10) 

i=occ 

The relation between the nuclear potentials and the meson fields will be derived in the 
next sub-section. In anticipation of these results, we set 

^^-3^>?(r-X;X), (11) 

and, for the simplest of the QMC models, 

^^^-2>9lpl{r-X-X). (12) 

The quantities p^(r — X; X) and p^(r — X; X) represent the nucleonic expectation values of 
the vector and scalar quark field operators for a nucleon with center-of-mass located at X. 
The variable r — X is the displacement between a quark and the center of its nucleon. For 
simplicity, we have assumed that the three quarks inside the nucleon have identical masses, 
and therefore identical density distributions. The calculation of these quark densities is 
described in the next sub-section. 

With these results, the meson field equations become: 

(-V^ + m^) u;(r) = ^gl j d'X p.(X)pl(r -X;X), (13) 
(_V2 + ml) a{r) = Sgj J d'X Ps{X)pl{r -X-X). (14) 

Thus the source term for uj{r) in Eq. (13) is a convolution of the vector density of the 
nucleon with the vector density of the quarks in the nucleon. Similarly, the source term for 



a{r) in Eq. (14) is a convolution of the scalar density of the nucleon with the scalar density 
of the quarks in the nucleon. This distinction, which was not present in the original work 
of Guichon [3], was first pointed out by Fleck et al. [4], who derived it from an analysis of a 
boost of the composite system. In Ref. [7] the convolution of densities has been replaced by 
an evaluation of the nuclear densities at the point r. This is equivalent to a local density 
approximation to the nuclear matter result, where the nucleon densities are constant over 
the volume of the nucleon. 

Strictly speaking, Eq. (13) is incomplete. The a;-meson couples to a quark baryon current 
of the simple form jj^ = ipq'^^ipq. However, the nucleonic matrix element of j^, which is the 
effective baryon current of the composite nucleon (■0) , will have both a 'ipj^ip term and an 
"anomalous" term dj.'t/ja^^tp. This is the well-known "tensor" coupling for vector mesons and 
photons. It is proportional to the anomalous magnetic moment of the nucleon, and emerges 
naturally in the present model. In the mean-field approximation there is a non-vanishing 
contribution from this term for the time-like part of the current — 0). We give a more 
detailed discussion of this point in Appendix A. Because the tensor coupling of the a;-meson 
is known to be small, we choose to ignore its influence in the computations performed here. 

The essential feature of the present sub-section is that the self-consistent solution of Eqs. 
(1), (13) and (14) give the energy and wave function of the nucleus. These equations are 
similar to the ones of QHD [6], but with the essential difference being the dependence on 
the quark vector and scalar densities. 



B. Quark Wave Functions and Densities 



Computing the quark wave functions and the resulting nucleonic vector potential \v{X) 
and effective mass M*{X) starts with using a fairly general representation of the field 
equation of a confined quark. 



d 

Z7°— -i-y ■Vr + mg + Vcon{r - X) 



ijUt,r-X) = 0. 



(15) 
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Here -X" is the position of the center of the free nucleon, niq is the bare quark mass, and 
the confining potential is defined as Vcon{f — -X"). This potential depends on the distance 
between the quark and the center of its nucleon. This central confining potential can be 
understood to represent the MIT bag, when complemented by the bag boundary conditions. 
In relativistic potential models one can take Vcon{v) — Cr^, as did Tegen et al. [9], who 
consider n = 2 and n = 3. Vcon can also be obtained as a Hartree approximation to a more 
realistic treatment of confinement based on two or three quark interactions. The superscript 
/ is meant to denote that the nucleon is free. 

Next we consider how Eq. (15) is modified when the nucleon is bound. Each quark feels 
a vector V^{r) and scalar V/(r) potential. If the nucleon is centered at a position X, then 
one expects that the quark wave functions will depend on both r and X. Thus we obtain 



^,(t,r-X,r) = 0. (16) 



This dependence can be fairly complicated. The ui and a fields and resulting quark potentials 
are functions of r, so that such fields are not functions of r — X, and are not central in the 
frame of the nucleon. Hence the ground state of a bound nucleon will not be spherically 
symmetric in general. 

To handle the difficulties that this entails, we develop a new approximation scheme based 
on the notion that 

Vv'M ij,{t, r-X,r)^ U,,s{X) 0,(i, r-X;X), (17) 

where C/y,s(X) is some suitable averaging of V^^g{r) over the volume of the nucleon. The 
notation (f)q{t, r — X; X) is meant to denote that the quark wave function depends on the 
coordinate variable r — X, and is only an implicit function of X through Vy^s{X). This 
notion is plausible because ipq{t,r — X ,r) vanishes if r — X is much bigger than the mean 
radius of a nucleon, and V7(r*) and V^{r) are not expected to vary much over this region. 
The variation of the nuclear fields is governed by two scales — the nuclear radius, and the 
nuclear skin thickness (distance at the surface over which the density decreases from 90% to 



10% of its maximum value) of about 2.5 fm. Either distance is much larger than the radius 
of the nucleon, so that replacing and by some suitable average over the volume of the 
nucleon seems reasonable, at least at the start. 

The procedure (but not the theory) of Ref. [7] is to evaluate the external fields at the 
nucleon center and neglect the variation of these fields within the volume of a single nucleon. 
This means that their calculations replace the approximate sign in Eq. (17) by an equality 

withc/„,,(x) = K%W- 

Here we intend to develop a more complete formalism, based on the plausibility of Eq. 
(17). We define an "average" potential Uy^s{^) and a residual interaction AUy^s{f, ^) such 
that 

K%ir)^U,,s{X)+AU,,,{r,X). (18) 

This equation is clearly a tautology. However, we may choose the potentials C/y,s(X) so as 
to make the first-order perturbation theory evaluation of AUy^s{f, X) vanish. This suggests 
that we can develop a convergent perturbative treatment of AUy^g- 

The feature that U^^g depends only on X allows us to simphfy the solution of the Dirac 
equation (16). We take 

^,{t, r-X,r)=J2cr, ^-iK(x)+u.(x)]t ^^^^ -X;X), (19) 

n 

where the quark energies e^{X) and wave functions (f)n{r — X; X) of quarks with quantum 
numbers n are determined by solving the equation 

l-zj -Vr + rrig- Us{X) + K„,(r - X)] 0„(r -X;X)= e*^{X) 7Vn(r -X;X). (20) 

The quantity c„ is to be determined by perturbative calculations with the residual interaction 
A[/,,,(r,X). 

Comparing Eq. (20) with the equation for the free nucleon (15), it is clear that the quark 
wave function in the medium is of the same form as the free quark wave function 0'^, but 
with the bare quark mass niq replaced by an effective mass 
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ml{X) = mg-U,{X). (21) 

Thus the quark energies and densities (or wave functions) are imphcitly functions of X only 
via m*{X). This is properly expressed by the notation e*{m*{X)) and pl g{r — X; m*(X)). 
We will continue to use the simpler notation e*{X) and pl g{r — X] X) with this understand- 
ing. We emphasize that the vector interactions have no effect on the nucleon properties in 
the medium other than an overall phase in the wave function, which results in a shift in the 
nucleon energies. 

As a first step let us assume that the nucleon consists of three quarks each in a state with 
K — —1. This state is treated in perturbation theory, so that state is simply 0o, which we 
shall simply denote as 0, i.e. 4> = 4'o- Similarly e = eo. (In this case c„ = c!^^ = 5n,o-) These 
ground state quark wave functions then determine the quark vector and scalar densities 

p?(f;X) = 0t(^;X)0(f;X), (22) 
p1{r;X)=^{f;X)<p{f;X), (23) 

subject to the normalization condition 

j d^f pl{f-X)^l. (24) 

We also introduce for convenience the definition 

jd'fpl{f-X)^S{X). (25) 

C/t, s(X) is chosen according to the criterion that the first-order perturbation theory 
evaluation of At/i,^s(r, X) vanishes. This is achieved with the definitions 

U,{X) = I pl{r - X- X)V,\r), (26) 
Us{X) = J d'r p%r - X; X)V,%r) / S{X). (27) 

Explicitly then, the first-order corrections to the potential Uv{X) and the quark energy 
e*(X) are 
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C/W(X) = I ct>^{r - X- X) AU,{r, X) ^{r - X; X) = 0, (28) 
^jj^^i ^ j _ j^. AUs{r, X) 0(r - X; X) = 0. (29) 

The quantity AU^^s does cause a first-order change in the quark wave function and a second 
order change in the potentials. This is discussed more fully in Appendix B. 

The net result is that Eqs. (26) and (27) determine the mean fields that act on the quarks. 
Our mean fields are computed by taking the average over the relevant quark densities, 
whereas in Ref. [7] these fields are evaluated at the center of the nucleon. This difference 
has a modest effect on the properties of finite nuclei, as discussed in the results described in 
the next section. 

We may now obtain the nuclear vector and scalar potentials to be used in Eq. (1). The 
term [/„ is present for each quark. Hence the nuclear vector potential is three times the shift 
in the quark energy caused by Uv{X), i.e. 

K(X) =3^7,(X). (30) 

The form of the effective mass M*(X) depends on whether one uses the relativistic 
potential models or the MIT bag model. In the relativistic potential models we expect the 
mass of the nucleon to be just the sum of the energies of the three quarks. Hence, in the 
absence of any correction for center-of-mass effects, 

M*(X) = 3e*(m*(X)), (31) 

where e*(m*) is the energy eigenvalue of Eq. (20). It is convenient to introduce a nuclear 
scalar potential Vs{X) to be used in Eq. (1), with 

VsiX) = M- M*{X) = 3 (e(m,) - e*(m*)) . (32) 

This facihtates a comparison of QMC models with QHD [6], and also allows us to use the 
experimental proton and neutron masses in Eq. (1) with the same scalar potential. 

In the MIT bag model with quarks of effective mass m*, the ground state solution to Eq. 
(20) is [2] 
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0(r) = 



joixqr/R) 
la- ■ fpqji{xgr/R) 



Xq 



(47r)V2' 



(33a) 



where 



— 



xi + {Rm 



1/2 



(O! - i?m!)/(0! + i?m!) 



1/2 



(33b) 
(33c) 
(33d) 
(33e) 



and Xq is the quark spinor. The eigenvalue Xq is determined by satisfying the hnear boundary 
condition jQ{xq) — (3qji{xq) at the bag surface. 
We take the nucleon mass to be 



M*{X) 



^ A fj 

' + ^ttR^B. 
R 



(34) 



Here B is the bag constant and 2; is a free parameter that is supposed to account for zero- 
point motion, vacuum corrections, etc. Our definition is different than Fleck et al. [4] and 
the initial work of Saito and Thomas [2], who include a correction for spurious center-of- 
mass motion. The parameters B and z are determined by the free nucleon mass M for a 
given bag radius Rq, and are given in Table I. A possible medium-dependence of B and z 
is considered in section IV. 

Following the Born-Oppenheimer approximation, the bag is assumed to respond instan- 
taneously to changes in the nuclear environment. Hence the equilibrium condition, 

dM* 

^R R=R* 

also applies in the medium [2]. In practice, for the models considered in this paper, R* at 
nuclear matter densities is not very different from the radius Rq of a free nucleon bag. Thus 
the dominant contribution to Vs{X) is given by the single particle energies, consistent with 
Eq. (32). 



0, 



(35) 
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The desired results of this section are the vector and scalar potentials that act on a 
nucleon, contained in Eqs. (30) and (32). These in turn depend on the mean vector and 
scalar quark potentials, given in Eqs. (26) and (27). The next step is to obtain more specific 
forms for these potentials for use in our calculations. This is done in the next section. 

C. Specific Evaluation of Nucleonic Potentials 

We wish to obtain the external vector and scalar mean-fields which act on quarks in 
a nucleon centered at X. We shall first consider the vector potential K(^)- The formal 
solution to Eq. (13) is 

//. p-mv\r-r'\ 
d'X' p„(X') J dy i^^^r^^ Pl{r' - X'; X'). (36) 

This is easily evaluated in momentum space. Using the relation 

47r|r-r'| " J (2^ q^ + ml' ^^^^ 
and changing variables to f = r' — X', we can rewrite Eq. (36) as 

V.\r) = 3{9tr J ^ / d'X' e--^V.(X) / d^f e-->2(f ; X'). (38) 

The last integral above suggests that one may define a vector form factor for bound nucleons: 

v{q;X)^ Jd're-^'^-'plir^X). (39) 

Note that the normalization condition (24) ensures that v{0;X) = 1. The use of the 
definition (39) leads to the result 

K'ir) = ^{9ir I ^ / d'X' e---'p.(X>(g; X'). (40) 

The external vector potential felt by a nucleon centered at X is given according to Eqs. 
(26) and (30) by the convolution 



K(X) ^sjd'r plir - X; X)K\r) 

d^q e*"-^ 
(27r)^ + ml 
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This quantity is easily recognized as the mean-field vector potential for nucleons interacting 
via meson exchange with (position-dependent) form factors. For the vector meson, the form 
factors are all that distinguishes the present model from the QHD model [6,10]. 
The quark scalar potential follows analogously from Eq. (40). We find 

K%r) = 3(,f)^ / ^ / d^X' 'p,(XO.(g; X% (42) 

where the scalar form factor s{q; X) for bound nucleons is defined by the equation 

s{q;X)^ j Sfe-^^-'pl{r-X). (43) 

Unlike the quark vector form factor, there is no normalization constraint on s(0; JC), which 
is just the quantity S{X) of Eq. (25). 

In analogy to Eq. (41), the mean quark scalar potential of Eq. (27) is given by the 
convolution 

U.(X) = j i\ pl(r - X; X)V^(t) I S(X) 

This is the quantity used in the Dirac equation (20) that determines the quark energies e*, 
which in turn determines the external nuclear scalar potential Vs{X). 

The total energy £ of Eq. (2) can now be given in terms of the quantities we have 
just determined. Using Eqs. (3), (13), (5), (40) and (41), we find for the vector meson 
contribution to the total energy 

E,^-\j d^X p,{X)V,{X). (45) 

This is the same as in QHD, and follows from the linear dependence of Vv{X) on the omega 
field oj{r). The effect is to remove half the vector potential energy contribution contained 
in the nucleon energies (sec Eq. (1)). 

For the scalar meson contribution, we use Eqs. (4), (14), (6), (42) and (44) to find 

Es^\j d^X p,{X) 3Us{X)S{X). (46) 
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Because Vs{X) ^ 3Us{X)S{X), this contribution does not remove half the scalar potential 
energy contribution contained in the nucleon energies. However, this relation is approx- 
imately correct, and becomes exact in the limit of large quark mass rUg. In this limit, 
S{X) 1, and Vs{X) becomes a hnear function of o-{r). One therefore expects to recover 
the results of the QHD model. 

Finally, we limit our considerations to spherical nuclei, so that both the quark and 
nucleon densities are spherically symmetric. The equations we need are then: 

V,{X) = / j,(qX)v{q; X) / dX' X'' 3,{qX)p,{X)v{q- X), (47) 

TT Jo + Jo 

Us{X)^-§-— dq^^jo{qX)s{q;X) dX X'' jo{qX')p,{X')s{q;X'), (48) 
OD ( A J TT Jo q^ + rrig Jo 

poo 

v{q; X) = AiT dr 3o{qr)pl{r; X), (49) 
Jo 

poo 

s{q;X)^An dr jo{qr)pj{r; X), (50) 
Jo 

together with the Dirac equations for the quarks (20) and for the nucleons (1). We have 
introduced the meson- nucleon coupling constants gy = Sg^ and gs = Sg^. 



D. Nuclear Matter Limit 

In the spirit of the usual approach to the QHD model [6] , we use the equilibrium proper- 
ties of isospin symmetric nuclear matter to fix the parameters of the model. If we consider 
the limit in which the nucleus is very large, the nucleon densities Pv{X) and Ps{X) can be 
treated as constants {pb and ps), as can the vector and scalar fields. Hence U-t, — and 
Us — V/, so that m* — mq — V^. Our formalism will yield the nuclear matter results of Refs. 
[2,4], except for the correction due to the spurious motion of the center-of-mass. According 
to Ref. [7], this correction should not be included. 

The composite nucleon obeys the Dirac equation with an effective mass M*, determined 
by Eq. (34), and a vector potential K = 3Vy^. Hence the nucleons have energy eigenvalues 

E{k) = {k^ + M*2)^/2 ^ ^yq_ (5;^) 
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The total energy £ is the sum of the energies of all nucleons below the Fermi energy together 
with the contribution from the meson fields, 

^-^r^s [(^^ + M*Y^ + 3K^] + - \mlu\ (52) 

Minimizing E with respect to uj and a yields the meson field equations (cf. Eqs. (13) and 

(14)) 

-9>=^-^ PB, (53) 

V! ^gla=^ p,S, (54) 

with 

Jo (27r)3 (A;2 + M*2)i/2' ^^^^ 

The last result for the quantity 5" is an MIT bag model result [2] . The nuclear matter results 
are therefore determined by the solution to the transcendental equation (54). In terms of 
the redefined constants — ?>g1 and gg — 3g'|, we can write 

£: = 4 ^ {k^ + M*^)V2 + A_pIs^ + A_pI. (57) 
Jo (27r)3 ^ ' 2mY 2ml 

The expressions (53), (54), and (57) depend only on the ratio of coupling constants to 

masses for both the omega and sigma mesons. Following Serot and Walecka [6] , we introduce 

the two dimensionless constants 

Cl^'^t (58) 



as free parameters. These are adjusted to reproduce the equilibrium properties of nuclear 
matter. We note that the present model reduces to the Walecka model in the hmit that the 
quark mass becomes very large while keeping the bag radius Rq fixed. This is a useful check 
on our numerical codes. 
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III. RESULTS 



In this section, we will use the MIT bag model to describe the nucleon. The two free 
nuclear matter parameters Cl and are adjusted to reproduce a binding energy per nucleon 
ol £ / pB — M — —15.75 MeV at a Fermi momentum kp — 1-30 fm~^. This corresponds to a 
saturation density of ps — 0.148 fm~^, which is the same as that used by Serot and Walecka 
[6]. A somewhat larger value of pb = 0.17 fm~^ (kp = 1.36 fm"''') was used by Saito and 
Thomas [2] . Our choice follows from the results of our finite nucleus calculations. With the 
larger saturation density we were unable to simultaneously fit the rms charge radii of light 
and heavy nuclei, whereas with the smaller value we could. 

Our nuclear matter results are given as models A1-A3 in Table II, corresponding to bag 
radii Rq of 0.6 fm, 0.8 fm, and 1.0 fm. The quark mass rUg has been set to 5 MeV for both u 
and d quarks. The bag radius R* decreases slightly in-medium, although it should be noted 
that the rms radius actually increases by about the same margin. In comparison with the 
results of the QHD model [6] , the effective mass is somewhat larger and the compressibihty 
somewhat smaller for QMC models. Models Bl and B2 allow for a variable bag constant, 
and model C uses relativistic oscillator wave functions instead of bag wave functions. These 
models are discussed in detail in the next section. 

For finite nuclei, the potentials Vy{X) and Vs{X) are used in the Dirac equation (1) 
for the nucleon to solve for the nuclear vector and scalar densities self-consistently. Our 
numerical procedure is to follow the standard iterative algorithm used in previous work 
[10]. Although the in-medium form factors v{q;X) and s{q;X) depend on m*{X), it is 
unnecessary to compute them at each value of X and again at each iteration. Instead, we 
pre-calculate these quantities for quark effective masses corresponding to Ug ranging from 
to 250 MeV in steps of 1 MeV. Intermediate values of m* are then obtained by interpolation. 

At this point we extend our discussion to include the p-meson and the photon (Coulomb 
interaction). The vector potentials Vp{X) and Vy{X) will have the same form as Eq. (47), 
with the replacement of the isoscalar vector density Pv{X') by the corresponding isovector 
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density for the p-meson, and the proton density for the photon. Because we have given the 
u and d quarks the same mass, the quark vector form factor v{q;X) is the same for the p 
and photon as for the lu. 

The p-nucleon couphng constant may be taken from experiment, or from the requirement 
that we reproduce the symmetry energy in nuclear matter [6]. Alternatively, it follows 
naturally within this model to assume the universal coupling = Qy, in which case gp = g^/S 
(since the p couples to the isospin of the quarks) . Because the a;-nucleon coupling constant 
g^ is much smaller here than in QHD, all three approaches give a value of gp that is roughly 
the same. We choose to fix gp = 2.63, which is consistent with the experimental value. (Note 
that the definition of gp used in refs. [2,6] is twice as large as ours.) 

As a check on our numerical codes, we were able to reproduce exactly the results of the 
Walecka model [6,10] in the limit that the quark mass becomes very large while holding the 
bag radius fixed at R* — Rq. 

Binding energies and rms charge radii for ^^O, ^'^Ca, ^°Zr, and ^°^Pb are shown in Table 
III. The scalar meson mass has been adjusted to fit the experimental rms charge radius 
of ^°Ca. The binding energy per nucleon for all QMC models shows some improvement 
over QHD. We note the strong positive correlation between binding energy per nucleon and 
M*/M at nuclear matter saturation density (see Table II). Charge radii are in reasonably 
good agreement with experiment for all nuclei. The charge density of ^^O for model Al is 
compared with experiment in Fig. 1. Other QMC models and QHD give similar results. 

The improvement over QHD in the binding energies comes at the expense of a reduction 
in the splitting of spin-orbit pairs, as shown in Fig. 2. This is easily understood if we recall 
that the spin-orbit splitting depends on the sum of the vector and scalar potentials [6] . One 
can therefore expect a strong negative correlation between spin-orbit splitting and M* /M, 
as seen in Fig. 2. We also note the insufficient binding of deeply bound states. 

Thus one cannot fit both the binding energy and single-particle energies within the 
framework of the Hartrcc approximation for this class of models. This points to the need 
for additional contributions to spin-orbit effects, such as tensor correlations. 
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IV. EXTENSIONS OF THE MODEL 



A. Bag Model 

One of the possible extensions of the QMC model is to allow the bag constants B and 
z to acquire a medium-dependence. The parameter z is supposed to incorporate ccntcr-of- 
mass corrections and the Casimir effect, both of which should depend on the quark effective 
mass m*. Initially we considered making z an arbitrary linear function of Us- However, 
for reasonable parameter ranges changing z had little effect, and tended to make the model 
worse instead of better. 

Alternatively, we can allow the bag constant B to vary in the medium. Since quarks 
presumably become deconfined at high enough densities, this would imply B ^ 0. This 
suggests a model for which B decreases with increasing density. A similar consideration has 
been made in a recent preprint by Jin and Jennings [8], who made B a function of M*. 
These authors also showed that for a particular functional form of B{M*), M* becomes a 
linear function of the scalar field a. Hence they obtained exactly the QHD model results for 
nuclear matter. 

For our purposes it is more convenient to assume that i? is a linear function of Us{X), 
such that 



B*{X) = B 



1 — as 



(60) 



M 

with ctfi an arbitrary parameter. Although B* does not go to zero at high density, we are 
restricting our consideration here to densities near that of nuclear matter. Typically in 
nuclear matter, Us ~ 130 — 200 McV for the models considered here. 

The additional dependence of M*{X) on Us{X), and therefore on cr(r), modifies the 
source term of expression (14) arising from Eqs. (8) and (12). Instead, we have 

as 



(61) 
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which follows from Eqs. (20), (27), and (34). Thus the effect of this additional term can be 
incorporated by making the change 

s{q;X)^s{q;X) 

to the results of section II C, and a similar change to the corresponding nuclear matter 
quantity S appearing in Eq. (54), viz. 

(63) 

The effect of B* is to reduce M* and increase the bag radius R* in-medium. This is 
reflected in the results shown in Table II for models labeled Bl and B2, corresponding to 

= 1 and = 2, respectively, and Rq — 0.6 fm. With increasing the results and 
parameters of the QMC models move closer to those of the QHD model [6] . 

For finite nuclei, the results also move towards those of the QHD model [6,10]. There is 
an improvement in spin-orbit splitting and in the binding energy of deeply bound states, as 
shown in Fig. 2. Model B2, with the largest change in the bag constant B*, does particularly 
well for the single particle energies and the charge density of ^^O, shown in Fig. 1. As noted 
in the previous section, an improvement in spin-orbit splitting comes at the expense of a 
modest reduction in total energy (Table 111). 

There are, however, differences between this model and QHD. All QMC models will have 
electromagnetic form factors that change in the medium. Indeed, the charge density, or at 
least the quark-core contribution to the charge density, is just the quark vector form factor 
v{q;X) (see Appendix A). This change arises principally from two sources. First there is 
a relative shift in the importance of the upper and lower Dirac wave function components 
for quarks due to the change in effective mass m*. Secondly, there is a change due to the 
change in bag radius R*. For models Bl and B2, the change in R* is quite large. We will 
explore the implications of medium-dependent electromagnetic form factors in a subsequent 
publication. 



1 + IttR* 



'B 



Ob 



3MS{X] 



(62) 



S ^ S 



1 + 



Ur*^b 



SMS. 
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B. Relativistic Oscillator Model 



As an alternative to the bag model, we consider a relativistic oscillator model with a 
(scalar) confining potential Vcon{r) = Cr'^, with C = 830 MeV/fm^. This model was used by 
Tegen et al. [9] to investigate the electromagnetic properties of the quark-core in nucleons. 
We refer the reader to this paper for further details. 

Our nuclear matter results given in Table II are quite similar to those of bag model Bl. 
In particular, the increase in the rms radius of the nucleon is of the same order as that for 
model Bl, which in the latter case is due predominantly to the decrease in the bag constant 
B. The results for the binding energies and rms charge radii are given in Table III. These 
results are again similar to those of model Bl, as is the spin-orbit splitting shown in Fig. 2. 

The relativistic oscillator model therefore offers an alternative to the bag models, and 
has one less free parameter than bag models with a medium-dependent bag constant B. Of 
course, we could also consider making the oscillator constant C medium-dependent. Another 
advantage of the oscillator models is that center-of-mass effects (which were considered by 
Tegen et al. [9]) are easier to handle than in bag models. 

V. SUMMARY AND OUTLOOK 

In this paper we have extended previous work on Quark-Meson Coupling models [2-4] 
in nuclear matter to finite nuclei. A generalized and systematic method for corrections to 
the model arising from the spatial nonuniformity of the meson mean fields over the volume 
of the nucleon is presented. Our procedure is to define the nucleonic potentials in terms of 
averages over the volume of the nucleon, as in Eqs. (26) and (27). For this definition, the 
polarization corrections are very small. 

The most basic versions of this model, based on either the MIT bag model [5] with a 
fixed bag constant or a relativistic oscillator model [9], have a small compressibility Ky^ 
compared with the QHD model [6], indicating a softer equation of state. While this has 



20 



experimental support, the drawback for these models is that this decrease in the compress- 
ibihty is accompanied by an increase in the effective mass ratio M*/M, which significantly 
reduces the spin-orbit component of the nuclear mean field and results in insufficient binding 
for deeply bound states. 

One possibility we have explored is to allow the bag constant B to decrease in the 
medium, reducing the effective mass M* and restoring much of the phenomenological success 
of the QHD model for finite nuclei [6,10]. The accompanying increase in the bag radius has 
implications for a number of properties of the nuclei, including electromagnetic effects such 
as charge densities, form factors, and the EMC effect [1]. 

In this work, we have ignored any effect of a decrease in the quark effective mass on 
the masses of the virtual mesons. Because the parameters we use depend only on the ratio 
of coupling constants to meson masses, a change in the meson mass will only result in a 
rescaling of the coupling constant, without affecting the nuclear matter results. However, 
there may be significant effects for finite nuclei, particularly at the nuclear surface, and we 
will explore these in the future. 

We have completely neglected the center-of-mass effects for quarks in the nucleon. It is 
important to handle this problem. In the light front formalism, the nucleon wave function 
is defined so that the intrinsic wave function is independent of the total momentum of the 
nucleon. Therefore, constructing a light front version of this model is a useful goal. 

There are other possible extensions of the quark coupling model. Aspects of chiral 
symmetry have been ignored here. It should be possible to construct a chiral version by 
using the cloudy bag model of the nucleon [12] in which the three-quark bag is surrounded 
by an evanescent cloud of pions. Yet another aspect is that correlations between nucleons 
are known to be important. Hence the Quark- Meson Coupling model should be extended 
beyond the mean field approximation. 

The construction of a chiral version of the quark meson coupling model, including physics 
beyond the mean-field approximation, and in which the nucleon center-of-mass degree of 
freedom is treated properly, is a useful and important goal. One would have a model 
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consistent with conventional nuclear physics, which embodies a reasonable treatment of the 
quark degrees of freedom. 
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APPENDIX A: ANOMALOUS VECTOR MESON COUPLING 



Consider the coupling of a vector meson to a quark current of the form — Qi'ipq'y^ipg, 
governed by some generahzed quark-meson couphng constant c/f to a quark of flavor i. We 
are interested in computing the effective "charge" and "anomalous" form factors governing 
the coupling of this meson to the composite nucleon. In the momentum representation, the 
matrix elements of the quark current between nucleon states \N{p) > may be written as 



N{p') 



i=l,3 



N{p))=gU{p')A^^{p'-p)U{p), 



(Al) 



where g is the meson-nucleon coupling constant, U{p) denotes a nucleon spinor of four- 
momentum p, and A'* is a vertex function of the form 



(A2) 



For quarks of equal mass, we can relate the coupling of the omega and rho mesons to the 
corresponding isoscalar and isovector parts of the familiar electromagnetic nuclear current. 
It is convenient to introduce the Sach's form factors 



GE{ql) = F,{ql)-r^F,{ql), 
GM{ql)^F^{ql) + F,{ql), 



(A3) 
(A4) 



with T) = — g^/4M^, and qf^ the square of the four-momentum transfer. For individual quark 
wave functions of the form 



0(r) = 



ig[r 
\ o- • rf{r) ] 



Xq 



(47r)V2' 



(A5) 



the expressions for the proton form factors are easily derived to leading order in 



GW)- / drr^ 3^{qr)\g\r) + f{r) 

Jo L 

Gniq^) = / drr^ ji{qr)g{r)f{r). 



(A6) 

(A7) 
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Note that G^O) = | {Gf{0) + g'^e\^)) = 1 and ^1(0) = | (^^^(O) - G^e\'^)) = 0, so that 
G^E^\o) = -Fi'^'^''(0) = 1 for both the isoscalar and isovcctor form factors. The magnetic 
moment of the proton due to the quark core is given by 

P Op roo 

/^P=^GUO) = --y^ drT'g{T)f{rl (A8) 

and the standard quark model ratio /x„ = — f/^p gives the magnetic moment of the neutron. 
The form factor G^^\q.) is just the quark vector form factor v(q) defined in Eq. (39). We ex- 
pect Qoj — for the (isoscalar) a;-nucleon coupling constant, and gp = gj for the (isovector) 
p-nucleon coupling constant. 

Experimentally, Fi'^\o) = -.1203 and F2^^^(0) = 3.706 for the electromagnetic form 
factors. Although this includes contributions besides the quark core, such as the pionic 
cloud, we expect the "anomalous" coupling of the vector mesons to the quark core to be 
of the same order as ^2°'^^ . Because the isoscalar anomalous couphng term is so small, we 
have chosen to ignore it in the calculations presented here. 

For completeness, we give the appropriate modifications of our formalism when the 
anomalous tensor terms are included. With the introduction of a "tensor" potential Vt{X), 
the Dirac equation (1) for the nucleon becomes 

-Z7 • V + M*{X) + 7°K(X) + Z7% • XVt{X)] ^'^(X) = Ea'^MX). (A9) 

It is convenient to introduce the "tensor" density 

Pt(X) = Yl ^MXh'-f ■ XA{X). (AlO) 

i=occ 

Assuming the nucleon wave functions are defined analogously to Eq. (A5), we have priX) — 
2Y.iG,{X)F,{X). 

For spherical nuclei, the vector potential Vv{X) of Eq. (47) becomes 



' jo{qX')p,{X')F,{q-X') 



2 roo roo 

+ <in{<iX')PT{x-'^'^''''''^ 



2M 



(All) 
24 



and the potential Vt{X) is 



Vt{X) = g^^ r dq jiiqX)^^^^ H dX' X'^ 

71 Jo q^ + m^ 2M Jo 



j^{qX')p,{X')F,{q-X') 

(A12) 



Note that Vt{X) oc —dVy{X)/dX in the absence of a medium-dependence to the form 
factors. Finally, the contribution of the vector meson field to the total energy is modified 
from Eq. (45) to 

E, = -\j d^X [p,{X)V,{X) + pt{X)Vt{X)] , (A13) 

which again removes half of the contribution of the vector potential contained in the nucleon 
energies. 

APPENDIX B: POLARIZATION CORRECTIONS 

Here we consider the second order perturbative corrections to the vector potential Uy{X) 
and the quark energy e*{X) arising from the residual interaction AUy^s{f,X). We also 
examine the first order corrections to matrix elements. All of these corrections involve 
intermediate excited states of the nucleon, hence the title of this Appendix. 

We study the corrections involving only the vector interaction here, and show that po- 
larization effects are vanishingly small. The formalism and physics of the scalar terms is 
similar, so we expect the same negligible effects. 

It is useful to introduce the braket notation, so that 

<r\cf>{X)>=cf>{r-X,X), (Bl) 

in the coordinate space representation. With this notation, the residual interaction can be 
expressed as (see Eqs. (18) and (27)) 

AC/„(r,X) = K%r)- < (/>(X)|K^(r-)|0(X) >, (B2) 
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so that the first order perturbation theory correction to Uv{X) vanishes: 

5C/W(X)=<0(X)|AC/,(r,X)|0(X) >= 0. (B3) 

The second order correction is given from standard Rayleigh-Schroedinger perturbation 
theory as 

where N and n represent the nucleon and excited baryon states, respectively. The sum on 
i is over the three quarks in the baryon. The confinement of the quarks to the region near 
X enables a Taylor series expansion of AUy{ri, X) about the point Vi — X. Let us define 
a variable 

s(^) =ri-X, (B5) 

so that 

At/,(s«,X) ^ V,\X) - < <p{X)\V,%rMX) > +3« ■ VV,^{X) 

The first three terms of Eq. (B6) do not lead to an excitation of the nucleon and do not 
contribute to the sum on n of Eq. (B4). This is because the first two terms do not contain 
any dependence on s^^\ and the third term is a center-of-mass operator. The last term of 
Eq. (B6) can be simplified by regrouping the term s^^s^^^ into quadrupole and monopole 
terms: 



Q% + SaAs^'^-s^'^. (B7) 



The quadrupole term, when combined with the spherically symmetric nature of V^{X), 
leads to a quadrupole dependence on the nuclear coordinates. The average of this vanishes 
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for our spherical shell model. Using the monopole term of Eq. (B7) in Eqs. (B6) and (B4) 
leads to the result 

^\<N{X)\j:i-i3hs^'^ ■s^'^n{X)>\ ^ ,2 
SUi'\X) = E J /TTv? F(x\ V^t^.(X) . (B8) 

n^N hN[X)-hn[X) ^ ^ 

This is an attractive second order contribution. 

We shall proceed to obtain a numerical estimate of the size of the correction given by 
Eq. (B8). First we shall assume that the isoscalar monopole excitation strength is taken up 
by a single state — the Roper resonance. This allows us to replace the different values of 
E]sf{X) — En{X) by a single one, AEr a; —500 MeV. The sum over n can then be performed 
using closure, so that 

SUi'HX) ^ ^^^^^ < N{X)\ < / > +A^- < >)' \N{X) >, (B9) 

in which the Jacobi coordinates p = ■^{s'^^^ — s^^^) and A = ■^{s'^^^ + s^^^ — 2s(^^) are used, 
and < p^(A^) > represents the nucleonic matrix element of (A^). Simple algebra, and 
the reasonable assumption that < A^ >=< >, which holds exactly in the hyperspherical 
approximation, leads to the expression 

SUi^\X) ^ ^^^^^ (< > - < >^) . (BIO) 

Dimensional analysis yields the result that < > — < >^oc r%, where tn is the root 
mean square radius of the nucleon. We proceed by using a simple Woods-Saxon form for 
t/,(X), 

U^{X) = Vol [l + e(^-^^)/"] . (Bll) 

Typically, Ra = tq A^^^, with tq pa 1.1 — 1.2 and a = 0.54 fm. 

The density is approximately constant near the center of the nucleus, so that 'V'^Uy{X) 
can be effective only near the nuclear surface. For example. Bit X — Ra 
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If one takes a and to each be of the order of the nucleon size, one sees that 5U!f^{RA) is 
proportional to (rjv/i?^)^, which is very small, a result that has been obtained without using 
a specific model of confinement. One may get a specific estimate by using the non-relativistic 
harmonic oscillator quark model of the nucleon. In this case, the closure approximation used 
above is true exactly, and < > — < >^= 1.5 rjf. Then, using typical values — 0.83 
fm. To — 1.15 fm, a = 0.54 fm, and Vq = 250 MeV, one obtains the result 

(i?^)^--J^MeV, (B13) 

which ranges from —0.5 MeV for ^^O to —0.1 MeV for ^°^Pb. These effects are truly small. 

We next consider the effects of p meson exchange. The equation that is analogous to Eq. 
(B6) is 

^ < 0(X)|l7(r)|0(X) > +r^) s« • VKf(X), (B14) 

in which the term 

s» is the dipolc operator £)*^*^. Now the term linear in s^^^ does 
contribute. We thus obtain the second order correction to the p exchange contribution to 
the vector potential: 

This term can be related to the experimentally measured dipole polarizability a, which is 
given by 



a 



|<iV(x)|E.=i,3^fW^)>' 



2e2 y I ^ V'^ \' \ ' . (B16) 



The use of the Wigner-Eckart theorem allows one to obtain the relation between SUf^ and 
the measured quantity a: 

(X) = -^ \VU,{Xt (B17) 

The use of the value a = 12 x 10~^ fm^ (as explained in the review [11]) and the Woods- 
Saxon form of the potential, along with the coupling constants discussed above, leads to the 
estimate SU^^\X — Ra) — —0.2 MeV, so that this effect is neghgible. 
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The previous examples show that the second order contributions to the vector potentials 
are very small. The use of the same formalism leads to the result that the second order 
contributions to the scalar potential are also negligible. 

Our final example concerns the computations of matrix elements. Consider some general 
observable O. The first order correction supplied by the vector potential is given by 

, ^ < N{X)\0\n(X) X n{X)\AU,\N(X) > 
To be specific, consider the case O = J2i s^^^ ■ s^*^ = s^- In this case we obtain 

or, using the approximations used to obtain the corrections to the energy, 

2 < > - < >2 
Hi^) = i AE 

In this case the contribution of the attractive scalar potential tends to cancel that of the 
repulsive vector potential. We account for this in a simple way by taking a net potential of 
depth 50 MeV. The net result is that 

< s2 > ^1/3 V ; 

which is again a negligible effect. 

The final result is that the dominant effects of the composite nature of the nucleon are 
contained in the definitions (26) and (27) of the mean quark vector and scalar potentials as 
expectation values of the quark vector and scalar potentials in the nucleon. 
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TABLES 

TABLE I. Nucleon bag model parameters corresponding to = 5 MeV and a free nucleon 

mass of M = 938.9 MeV. 

Ro (fm) ^ttR^B (MeV) z 

0.6 232.918 4.00367 

0.8 232.915 3.29545 

1.0 232.912 2.58724 



TABLE IL Model parameters corresponding to a nuclear matter binding energy per nucleon of 
—15.75 MeV at kp = 1-30 fm~^. Nuclear matter results are given for the compressibility, nuclear 
effective mass, and in-medium bag radius (rms radius for the oscillator model). 



Model 


i^o (fm) 






Ky^ (MeV) 


M*/M 


R* (fm) 


QHD 




357.62 


274.00 


547 


0.541 




Al 


0.6 


154.94 


117.46 


292 


0.775 


0.596 


A2 


0.8 


132.93 


98.41 


279 


0.803 


0.795 


A3 


1.0 


117.15 


84.57 


266 


0.823 


0.993 


Bl (ae = 1) 


0.6 


187.98 


145.50 


321 


0.734 


0.624 


B2 (ae = 2) 


0.6 


289.53 


226.11 


364 


0.614 


0.676 


C (Rel. Osc.) 


0.642^ 


190.09 


147.26 


326 


0.731 


0.689^ 



^For the relativistic oscillator we give the rms radius. 
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TABLE III. Binding energy per nucleon (in MeV), and rms charge radii (in fm) for several 



closed shell nuclei. We have taken 


m„ = 783 MeV, rup 


, = 770 MeV, and gp 


= 2.63. 


The scalar 


meson 


mass nis 


has been adjusted to fit the 


rms charge 


radius of 3.48 fm in '^^ 


'Ca. 




Model 


rUs 




"O 




^°Ca 




90 r? 


208pb 




(MeV) 


t/A 


^ 2^1/2 
< ^ >ch 


8/ A 


o 1/9 

< >1 


CIA 

£/A 


^ 2^1/2 
< ^ >ch 


8/ A 




QHD 


520 


-4.88 


2.7 A 


-6.30 


3.48*^ 


-7.07 


4.27 


-6.82 


5.48 


Al 


450 


-6.22 


2.75 


-7.46 


3.48^ 


-7.84 


4.32 


-7.59 


5.59 


A2 


415 


-o.yo 


z. ( 


7 


3.48^ 


7 70 




-7.49 


5.59 


A3 


385 


-5.72 


2.75 


-7.12 


3.48'' 


-7.56 


4.32 


-7.40 


5.59 


Bl 


460 


-5.93 


2.75 


-7.23 


3.48^ 


-7.68 


4.31 


-7.46 


5.57 


B2 


480 


-5.57 


2.75 


-6.94 


3.48^ 


-7.55 


4.29 


-7.34 


5.52 


C 


455 


-5.51 


2.75 


-6.94 


3.48'' 


-7.46 


4.31 


-7.30 


5.56 


Expt 




-7.98 


2.75 


-8.55 


3.48 


-8.71 


4.28 


-7.87 


5.50 


^Fit 
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FIGURES 

FIG. 1. The charge density of ^^O for bag models Al and B2, which both have Rq = 0.6 
fm. Model B2 has a medium-dependent bag constant. The solid line is the experimental param- 
eterization of H. de Vries, C. W. de Jager and C. de Vries, At. Data Nucl. Data Tables 36, 495 
(1987). 

FIG. 2. Proton single particle energies in ^^O for all models. 
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